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Abstract 

QCD sum rules are used to compute the first few moments of the mesonic quark momen- 
tum. Transverse, longitudinal and mixed transverse-longitudinal components are examined. 
The transverse size of the pion is shown to be dictated by the gluon condensate, even though 
the mass and the longitudinal distribution are dominated by the quark condensate. The 
implications of our results for color transparency physics and finite temperature QCD are 
discussed. 
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Mesonic wave functions have been the focus of much recent theoretical activity regarding 
exclusive processes [p], ^], and color transparency [Q. In the former case, the asymptotic 
formula of exclusive process can be expressed, via factorization, in terms of a hard scattering 
amplitude T H and the quark distribution function Q 2 ) || (£ is the light-cone momentum 
fraction of quarks and Q 2 is the typical external momentum). However, there have been 
many suggestions that such formulae are not valid at presently accessible energy ranges 
!| [5J. Furthermore, computations of form factors often employ a wave function with a 
spread of b (the transverse separation between quarks) |], |8|. Sterman and collaborators 
have computed the Sudakov effects §| for the pion and proton form factors[0, p7o|| . This 
involves "wave function" of the form e s<yb ^ , Q \ The Sudakov effect should dominate at high 
Q 2 , but it would be surprising if the 6-dependence resides only in the function 5*. 

The transverse size dependence plays also a key role in color transparency physics. A 
Gaussian dependence on the transverse momentum k\ (k±_ is canonically conjugate to b) 
implies that the effects of color transparency would never be observed, but a power fall off 



does allow such observations gj. Furthermore, QCD lattice calculations |1 1J now provide 
hadronic wave functions, and understanding these with analytic techniques would be useful. 

The use of QCD sum rules (QSR) has been a particularly useful example of such an 
analytic technique. QSR have long been known to provide reasonable estimates of hadronic 
properties |L2| and distribution functions needed as inputs to perturbative QCD (PQCD) 



calculations. In particular, Chernyak and Zhitnitsky [13, 14] have determined £ 2n moments 
of the distribution function </>(£). 

The present paper, concerns a new application of QCD sum rules(QSR): the determina- 
tion of the ^ 2l k 2 ^ moments of mesonic wave functions. We derive new Borel sum rules for 
these moments, which makes us possible to extract physical picture of the quark momentum 



distributions. (For an application of a finite energy sum rule, see ref. |15|.) The starting 
point is the Bethe-Salpeter amplitude for the pion 

if7rPa4>A z -P^ 2 )^ + z a xA z -V,z 2 )v = (0|d(2)7 a 7 5 e i!?j -2 a ' u(-z)\it + {$>)) ^ 

(j\ n _ <-> 

= E^m^iiz-DTum^ip)),, (i) 

n n - 

where <p n and x-t represent two different Bethe-Salpeter amplitudes, i D= iD + {iD)\ D = 
d — igB a ^-, n is the renormalization scale and is the pion decay constant. First let us 

introduce a light-like four vector y and a transverse vector t = (0, b, 0) which is perpendicular 
to the hadron momentum p = (p°,0,p z ). Then we make a decomposition z = y + t so that 
z 2 = —b 2 , and z ■ p = y ■ p. Since we have chosen z to have both longitudinal and transverse 
components, the matrix element in eq.(l) can be further expanded in terms of y and — b 2 . 
The result is that 

(0\d(0)y a la (iz- D)M0)\* + (p))» = iUit(y-p) n - l+ \-b 2 ^M n _^ (2) 

1=0 

after contracting with y a , and where the normalization is M 0i o = 1- Since y ■ z = 0, the con- 
traction of y a with eq. (1) enables us to extract (f) n on which we focus our attention in this pa- 
per. M n _ij are related to the matrix elements with n—l covariant derivatives in the longitudi- 
nal direction and / in the transverse direction. Both n and / are even numbers. The moments 
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M n -ij can be related to spin n + 1 and twist I + 2 matrix elements by identifying independent 
symmetric traceless matrices. For example, —if^Mo^P^ = ^i(0\d'y fl r y 5 <T al 3gG al3 u\ii( y p)) . 

The function </> n (z ■ p, z 2 )^ depends on coordinates; the relation to the momentum space 
wave function is given by 

MV ■ V, -b 2 ), = £ d£ f d 2 k ± e^e- ik - rb M£, k 2 ± U (3) 
while the light-cone quark distribution function is given by 

MZ)v = J d 2 k ± 6(fx 2 - k 2 ± )Mt, *?±V ( 4 ) 
Expanding the integrand of eq.([jj) as a power series in y ■ p and (b 2 )^ leads to 

n n - i=o 

where 

(C' l k l ± ), = £ / d 2 k ± C~ l k l ± M^ (6) 

= T(l/2 + Z/2)/(r(l/2)r(l + 1/2)) and n Q = T^pj. The desired moments are now 
identified by using eqs. @ and eq.(H) in the contraction of eq. (1) with y a and equating the 
coefficients of the terms in the expansion. The result is 

< C~ l k l ± >= -i-M^C-l)'/*. (7) 

n.W-L I 

We proceed by using QCD sum rules to extract M n -i t i from eq.(H) an d thereby obtain 
the moments from eq.(0). Consider the correlation function T n .o, 

T ni0 (y ■ q, -b 2 ,q 2 ) = i J d A xe^{0\T[d(x)y a lal5 (z- D) n u(x), u(0)y a lal5 d{0)}\0) . (8) 

The complete sum over intermediate states is approximated by the lowest mass pion pole 
term and a continuum contribution. Then the spectral density can be expressed as 

1 n 

-ImT nfi (y ■ q, -b 2 , q 2 ) = £ f 2 (y • g)^+ 2 (-6 2 )! (M n _ M % 2 - m 2 ) + c(n, l)8(q 2 - s n )) , (9) 
71 1=0 

where c(n, I) are coefficients chosen to match the perturbative part in the Operator Product 
Expansion (OPE) at Q 2 — > oo. We have carried out the OPE for T n for operators of 
dimension less than or equal to six: 

n 

T n , (yq,-b 2 ,q 2 ) = T.nCMy.pr-^ib^^l-CZM-e) 

1=0 

-KM-v 2 ) + <%)<~&)^ - cii(-i) ,/a <VSw> a £], (10) 
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where 

ci, = l B( rL^ + i4 + 1 ) ( !) ( ^) 

[1 + 9(n - 2 - /)(« - + »(i - 2)f^], 

Cl,, = ± B (^ + i,I + l) (I -l) 

[1 + 6(n - 2 - /)(n - ()^ + »(' - ^'jT^l' 

Ct, = +^o^(4n+ll) + ^ 2 g-g f -n(n-l) ) (11) 

with B(n,l) = r(n)r(/)/r(n + /) and #(n) = 1 for n > 0. The OPE for the p meson is as 
above except for the four quark condensate term for which (An + 11) — > (4n — 7). 

We may immediately examine the physical consequences of these expressions. The gluon 
condensates arising from the covariant derivative in eq.(|J) give the terms proportional to #'s 
in C\ i and G\ x . These terms together with the other gluon condensate become increasingly 
important as I increases. Conversely, one of the four quark condensates with hard gluon 



line (Fig. 3 of []13|1 ), which gives the dominant power correction for the I = sum rules, 
does not contribute at all for 1^0. This is because the internal gluon line in the relevant 
tree graph carries no transverse momentum. Thus the power correction is dominated by 
the gluon condensate and the quark condensate plays only a minor role. This is opposite to 
the usual sum rules for longitudinal moments (I = 0) and for the meson masses where the 
quark condensate term is essential. In the constituent quark picture of hadrons, the mass of 
the constituent quark is obtained through the dynamical breaking of chiral symmetry and 
the transverse wave function is determined by the confining force between the constituent 
quarks. This picture is consistent with our observations about the OPE results. 

To go beyond these qualitative aspects we follow the QSR standard procedure and, 
equate the phenomenological side eq.@ to the OPE side eq.fllQl) using the the standard 
Borel transformed dispersion relation on q 2 , - J dse~ s ^ M lmT n (s) = Borel trans [ReT n ], 
and identify the coefficients of the double expansion. The result is 

d_ 

d{i/M 2 



+^(2 - l^Uii-^T^iMr^i-G 2 } 



2 /" * 



+0(l-4)Cl +l /^)KMr 2 - 1 (-G 2 ) 

l 71 
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+0(2-/)^ +M (-l)'- O r T (M 2 )^ (v ^ g - g)2 . (12) 

I 2 )• 

The above sum rule reduces to the sum rule for f n [[12] for n = I — and to the 



longitudinal sum rule [13], [Tj] for / = 0. We will use the physical value f n = 133 MeV 
and derive each moment for different value of So- For the vacuum condensates, we use the 
standard values [ |T2]j : 



auu) 2 ~ 1.83 • 10~ 4 GeV 6 , (-G 2 ) ~ 1.2 • 10" 2 GeV 4 . (13) 

7T 

The outline of the method used to analyze the sum rule is the following: 

1. Find the Borel window, i.e., find M^ in such that the total power correction is less than 
30 % of the perturbative part. This will insure that the neglected higher dimensional 
operators cannot be large for M > M min . We then choose M^ ax arbitrarily say M^ ax = 
M Ln + °- 3 or 0.4GeV 2 . 

2. Find the value of So which minimizes the dependence on M 2 . Then the extracted 
moments will not depend too much on the exact choice of M^ ax . 

3. The physical quantity is then obtained from the average over the Borel window M 2 . 

Chernyak and Zhitnitsky used a similar procedure for the case of I = 0. 

The form of the Wilson coefficients dictates the characteristic behavior of the sum rule as 
a function of the variables n, I. Increasing n for fixed I generally enhances the coefficient 
of both the gluon and 4-quark condensate in the OPE. Then the Borel window (range of 
M 2 ) appears at higher values of M 2 and the value of So increases. Note also that the power 
correction becomes large as n increases, so that the sum rule is not reliable for n > 6. 
Increasing / for fixed n generally enhances the coefficient of the gluon condensate but 
will reduce that of the 4-quark condensate. In particular, there is no 4-quark condensate for 
/ > 2. This competing effect tends to keep the Borel window fixed or moves it to slightly 
higher values of M 2 . For I > 6, the sum rule is again sensitive to So and one cannot make 
reliable estimates. The value of So must be reduced compared to other moments to reduce 
the M 2 dependence. 

Table 1 shows the result of our analysis. The moments for I = correspond to the values 



obtained in ref. |13|, 14}] . The main sources of errors in Table 1 are the contribution from 
higher dimensional operators and the uncertainty in the exact value of So- The first issue 
is related to Mikhailov and Radyushkin's (MR) ||16[ criticism of ref. [13], [nj] . MR stress the 



importance of including non-local condensates [l7| . Although, there are some ambiguities 
associated with modeling different types of non-local condensates, we take the MR results 
as a guide to estimate the uncertainties of our calculation. The Gaussian model of MR leads 



to longitudinal moments about 50 % smaller than that of ref. [14]. We expect similar 
behavior to be true for sum rules with I ^ and estimate the uncertainty associated with 
unknown higher dimensional operators to be —50 % for all moments. A second source of 
errors is the lack of knowledge of So which strongly influences the knowledge of moments 
with / ^ 0. For example, increasing the value of So to its value for / = 0, causes a 50% 
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increase in the moments. This gives an upper limit for the I 7^ moments. So overall, the 
estimated errors for all the moments are ±50 %. 

In the conventional notation, where the separation between quarks is b instead of 2b, our 
n = 0, / = 2 moment implies an average transverse momentum of (300MeV) 2 . This is a 
characteristic hadronic scale and is consistent with a rough estimate (323MeV) 2 in ref . [|T4[ . 



We can use our results to study three features of the wave function; approximate factor- 
ization of the longitudinal and transverse directions; power law behavior of the transverse 
wave function; and the low temperature behavior of the transverse wave function. Let us 
discuss each one separately. 

Here the term "factorization" refers to the property that (£ n fc^} ~ (£ n )(^j_)- A recent 



calculation |]T8[ of color transparency effects in high energy pion-nucleus scattering assumed 
this factorization. The results shown in Table 1 are consistent with the factorization property. 
The numerical origin of factorization is natural and can be traced from the expressions 
provided here. Further work is in progress to demonstrate that this factorization holds when 
including the effects of non-local condensates for I ^ 0. 

The second point is the k± dependence of the transverse wave function. Once we assume 
factorization, the transverse momentum dependence can be examined using the relation 

0(£,^)=0(£M^)- (14) 

The knowledge of a finite number of moments is not sufficient to completely determine 
ip(k±); instead we define a quantity A such that A = \{k\) — (k\) 2 }/ (k\) 2 . We can then 
determine which of two popular forms, a Gaussian ipG{k±) = Aexp(— k\/m 2 ) and a power 
law ippL(k±) = A ^ k -2 ^2)2 , has a value of A consistent with the results of Table 1, which 
is that A ^ 8. One may use the ^'s to calculate A; it is first necessary to introduce a 
value of /i 2 (see eqs. (3) and (4)) which is expected to be of the order of 1 GeV 2 . For 
the Gaussian wave function ipG{k±), Ag ~ 1 — 0(exp(—fi 2 /m 2 )). The exponential term 
can be estimated by reproducing the value of (fc^) shown in Table 1. For /x 2 larger than 1 
GeV 2 , the 0(exp(—fi 2 /m 2 )) turns out to be always less than 0.1. For the power law, we can 
again determine [i 2 by reproducing the value of {kj_). Then Apl = 28 for /i 2 = 1 GeV 2 and 
Apl = 6 for fi 2 = 2 GeV 2 and always larger than 5 for larger values of /i 2 . So the sum rule 
values for A are more consistent with a power law type of transverse wave function even if 
we take into account the estimated errors. More generally, the concept of a large value of A 
is consistent with the notion of significant fluctuations of the momentum, a property that 
favors the possibility that color transparency would occur 

Another interesting question is the low temperature generalization of our moment sum 
rules. This generalization is achieved by applying the low temperature pion gas approxi- 



mation to the thermal expectation value of the hadronic correlation functions. [|19| , |20|| The 
use of soft pion theorems lead to the result that to lowest order in e = T 2 /6f 2 , the axial 
vector and vector correlation functions at finite temperature can be expressed in terms of a 
linear combination of vector and axial vector vacuum correlation functions with temperature 
dependent residues ||19|, such that the only effect at finite temperature is a renormalization 



of /J = f n (l — e/2). We find that the same is true for the transverse moments, i.e., 

M A t (T) = (1 - e)M A l (T = 0) + eM^T = 0) 
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M^(T) = (1 - e)M^(T = 0) + eM^T = 0), (15) 

where the superscripts V and A represent moments for the vector and axial vector currents. 
So the change at finite temperature can all be accommodated by a change in /„- with no 
change in the transverse moments. This suggests that to lowest order in T, the wave functions 
do not change when the temperature is increased from zero to a small value. 

In summary, we have derived sum rules for the transverse and longitudinal moments of 
the pion (and the p meson) wave function. These results can be used to understand the role 
of the quark wave functions in exclusive processes and at finite temperature. 

This work was supported in part by the US Department of Energy under grant DE-FG06- 
88ER40427. 
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Z = 


1 = 2 


Z = 4 


moments 


M 2 ■ 

mm 

(GeV 2 ) 


So 
(GeV 2 ) 


moments 
(GeV 2 ) 


M' 2 ■ 

mm 

(GeV 2 ) 


S 
(GeV 2 ) 


moments 
(GeV 4 ) 


M' 2 ■ 

mm 

(GeV 2 ) 


So 
(GeV 2 ) 


n = 


1 


0.6 


0.9 


0.36 


0.8 


0.5 


1.17 


1.2 


0.4 


n = 2 


0.4 


1.3 


1.9 


0.14 


1.3 


0.9 


0.30 


1.5 


0.5 


n = 4 


0.24 


1.9 


2.8 


0.09 


1.9 


1.4 


0.20 


2.0 


0.5 



Table 1: Mixed moments of the quark transverse/longitudinal momentum distribution. 
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